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O ■ ABSTRACT 

Context. The observed presence of pure CO2 ice in protostellar envelopes, revealed by a double 
peaked 15 fim band, is often attributed to thermally induced ice segregation. The temperature 
required for segregation is however unknown because of lack of quantitative experimental data 
and this has prevented the use of ice segregation as a temperature probe. In addition, quantitative 



segregation studies are needed to characterize diffusion in ices, which underpins all ice dynamics 
and ice chemistry. 



Aims. This study aims to quantify the segregation mechanism and barriers in different H20:C02 
f^*) ' and iLO:CO ice mixtures. 

Methods. The investigated ice mixtures cover a range of astrophysically relevant ice thicknesses 
and mixture ratios. The ices are deposited at 16-50 K under (ultra-)high vacuum conditions. 
Segregation is then monitored, at 40-70 K in the CO2 mixtures and at 23-27 K in the CO mix- 
tures, through infrared spectroscopy. The C0 2 and CO band shapes are distinctly different in pure 
and mixed ices and can thus be used to measure the fraction of segregated ice as a function of 
time. The segregation barrier is determined using rate equations and the segregation mechanism 
is investigated through Monte Carlo simulations. 

Results. Thin (8-37 ML) H2O ice mixtures, containing either CO2 or CO, segregate sequentially 
through surface processes, followed by an order of magnitude slower bulk diffusion. Thicker ices 
(> 100 ML) segregate through a bulk process, which is faster than even surface segregation in thin 
ices. The thick ices must therefore be either more porous or segregate through a different mech- 
anism, e.g. a phase transition, compared to the thin ices. The segregation dynamics of thin ices 
are reproduced qualitatively in Monte Carlo simulations of surface hopping and pair swapping. 
The experimentally determined surface-segregation rates follow the Ahrrenius law with a barrier 
of 1080 ± 190 K for H 2 0:C0 2 ice mixtures and 300 ± 100 K for H 2 0:CO mixtures. Though the 
barrier is constant with ice mixing ratio, the segregation rate increases with CO2 concentration. 
Conclusions. Dynamical ice processes can be quantified through a combination of experiments 
and different model techniques and they are not scale independent as previously assumed. The 
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derived segregation barrier for thin H 2 0:C02 ice mixtures is used to estimate the surface segrega- 
tion temperature during low-mass star formation to be 30+5 K. Both surface and bulk segregation 
is proposed to be a general feature of ice mixtures when the average bond strengths of the mixture 
constituents in pure ice exceeds the average bond strength in the ice mixture. 

Key words. Astrochemistry; Line: formation; Molecular processes; Methods: laboratory; 
Circumstellar matter; ISM: molecules 



1. Introduction 



Ices form in dark clouds through accretion of atoms and molecules onto cold (sub)micron-sized 
dust particles. The atoms and molecul es are subsequently hydrogenated and oxygenated on the 



grain surface to form ices such as H 2 (Me rrill et al.ll 1 976c iTielens & Ha gen 1982; 



2003 



seqi 



Watanabe et al 



Ioppolo et al.l l2008). This surface formation process is efficient enough that i n the densest 



star-f orming regions up to 90% of all molecules, except for H 2 exist in ice form (e.g. lBergin et al 



2002). Observations of dense clouds show that two of the most abundant ices, H 2 and C0 2 , 



(Bergin et al. 


2005; 


Pontoppidan 


2006; 



Sonnentrucker et al 



2008). These results predict a 



bi-layered ice composition, with a bottom layer dominated by a H 2 0:C0 2 mixture of ~5:1 and 
top layer consisting of CO-rich ice. Infrared observations of the C0 2 ice bending mode around 
15 yum towards dark clouds confirm this scenario; the C0 2 spectral band consists of two or more 
distinct componen ts consistent with laboratory ice spectra of H 2 0:C0 2 and CO:C0 2 ice mixtures 



Knez et al. 



2005). In pure CQ 2 ice spectra, the bending mode has a characteristic double peak 



S andford & Allamandola 



199ol e.g.). This double-peak is not observed towards star-forming re- 



gions before the protostellar collapse and the turn-on of the protostar. This agrees with current 
astrochemical models, where pure C0 2 can only form through thermal processing of previously 
mixed H 2 0:C0 2 and CO:C0 2 ices. 

Heating of ice mixtures in the laboratory results in both sequential desorptio n, starting with 



the most volatile molecules, an d in ice re-structuring, including ice segregation (e.g. lCollings et al 



2004 



Ehrenfreund et al 



1998). In H 2 0:C0 2 and CH30H:C0 2 ice mixtures, ice segregation is 
identified from the growth of pure C0 2 ice features, especially the characteristic 15 /urn double 
peak. Therm al desorption of CO from CQ 2 ice mixtures also results in the same pure C0 2 double- 
peak feature ( Ivan Broekhuizen et al Lood) . Ice heati ng is therefore inferred towards low and high- 



mass protostars wherever pure CQ 2 ice is obse rved dGerakines et alJll999 ; 



Nummelin et al 



2001 



Pontoppidan et al 



2008 



Zasowski et al. 



2009). Towards most of the protostars, low- and high- 



mass sources alike, the majority of the C0 2 ice is still present in the two phas es typical for dark 



cloud cores, suggesting a range of ice temperatures in the protostellar envelope dPontoppidan et al 



2008). The C0 2 spectra then contain information about the temperature structure of the protostellar 
surroundings. Maybe more importantly, the spectra contain information on the maximum temper- 
ature to which the icy grain mantles in the envelope have been exposed, since ice processing is 
irreversible. The fraction of pure C0 2 ice, with respect to the total C0 2 ice abundance towards a 
protostar, can therefore be used to constrain its thermal history, and its variability compared to cur- 
rent stellar luminosity measurements. Attempts to use this information quantitatively suffers from 
difficulties in determining the origin of the pure C0 2 ice, because of a lack of quantitative data 
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on H20:CC>2 ice segregation, i.e. the segregation temperature of H20:C02 mixtures at these time 
scales is unknown. 

Segregation is expected to occur whenever diffusion of molecules in the ice is possible and 
it is energetically favorable for molecules of the same kind to group together; for example H2O 
molecules form stronger hydrogen-bonds with each other than with CO and CO2. Thus, ice seg- 
regation studies provide information on diffusion barriers. These barriers are currently the most 
important unknowns in models of complex organic ice formation, where la rge molecules form 



through recombination of smaller molecules and radicals that diffuse in the ice (IGarrod et al. 



2008). 



Ice diffusion barriers also govern selective desorption from ice mixtures, which is responsible for 
much of the chemical structure in protostellar envelopes. While experiments on ice chemistry and 
desorption provide some information on diffusion, the effects of diffusion are difficult to disentan- 
gle because of the additional dependencies of the experimental results on desorption and reaction 
barriers. Segregation studies offer a comparatively 'clean' environment to quantify diffusion within. 



under high-vacuum conditions by multiple groups ( 


ihrenfreund et al, 


1998 


1999; 


Dartois et al. 


1999 


Palumbo & Baratta 


2000 


Bernstein et al. 


2005 


Oberg et al. 


200" 


). The results differ some- 



what between the different studies, but in general H20:C02 ice mixtures are observed to segregate 
between 60 and 75 K at laboratory time scales for mixtures between 9:1 and 1:1. Some of t he ob 
served differences are probably due to the different mixing ratios, since 



Ehrenfreund et al 



(1999) 



demonstrated that segregation in CH30H:C02 ice mixtures depends on the original ice mixture 
composition; at 60 K only ice mixtures with ten times more CO2 than CH3OH segregate, at 100 K 
CH30H:C02 1:3 ice mixtu r es seg regate efficiently as well. 

Recently iHodvss et alj ( 12008 ) investigated segregation in two 0.15 fim thick, H20:C02 ice 
mixtures (9: 1 and 4: 1 , respectively) under high vacuum conditions. They observed an onset in seg- 
regation at 60 K, which they interpreted as resulting from the known, slow H2O phase change from 
high density to low density amorphous ice between 38 and 68 K. The immediate segregation of 
an ice mixture deposited at 70 K shows, however, that segregation may also occur through dif- 
fusion. Consistent with previous inves tigations they found segregation to be more efficient in the 



4:1 compared to the 9:1 ice mixture. 



Hodvss et al 



(2008) also investigated the growth of segre- 



gated, pure CO2 in the 4: 1 ice mixture with time at three different temperatures and they observed 
a faster growth of the pure CO2 feature, and also an increased final segregated fraction, at higher 
temperatures. The ice deposition temperature was important above 50 K in this study, while lack 
of quantified uncertainties prohibited an estimate of its significance at lower temperatures. 

Building on these previous results, we aim to systematically investigate the segregation behav- 
ior of H20:C02 and H20:CO ice mixtures under astrophysically relevant conditions. The focus is 
on the CO2 mixtures because of their previous use as temperature tracers around protostars. The 
CO mixtures are mainly included to test whether the deduced segregation mechanism for H20:C02 
mixtures is molecule specific or if it can be applied more generally; The H20:CO segregation tem- 
perature and mechanism are however astro physically relevant as well since H20:CO ice is present 
in star-forming regions dChiar et alJll994l) . For both ice mixtures, the segregation rates are mea- 
sured under a range of experimental conditions, which are discussed in Section 2. The experiments 
are complimented by Monte Carlo simulations, introduced in Section 3, to test the theoretical out- 
come of different segregation mechanisms. The results of ultra-high- and high-vacuum experiments 
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on segregation are presented together with the deduced segregation barriers and the simulation re- 
sults in Section 4. Based on the simulations and the experimentally determined segregation-rate 
dependencies, the possible segregation mechanisms are discussed in Section 5 followed by astro- 
physical implications. 



2. Experiments 



Quantitative experiments are carried out with thin ices, < 40 monolayers (ML), under ultra-high 
vacuum (UHV) conditions (~10" 10 mbar) on the set-up CRYOPAD. Complimentary experiments 
on thick ices (>100 ML) are performed under hi gh vacuum (HV) co nditions (~10" 7 mbar) at a dif- 



ferent set-up. CRYOPAD is described in detail by 



Fuchs et al. 



(2006). The set-up is equipped with a 



Fourier Transform InfraRed (FTIR) spectrometer in reflection-absorption mode (called Reflection- 
Absorption InfraRed Spectroscopy or RAIRS), which covers 4000-8 00 cm" 1 with a typical s pectral 



Gerakines et al 



1995) and is 



resolution of 1 cm" 1 . The HV experiment is described in detail by 
equipped with a FTIR, set up in transmission mode, which spans 4000-600 cm" 1 and is run at 1 
cm" 1 resolution. 

In all experiments, the ices are built up diffusively on a gold surface (UHV) or a Csl window 
(HV) by introducing pre-mixed gases into the chamber at the chosen deposition temperature. The 
gas mixtures are prepared from 13 CC<2 and 13 CO (Indugas, 98-99% isotopic purity) and deionized 
H2O, which is further purified by several freeze-pump-thaw cycles. Because all experiments are 
with 13 CO and 13 CC>2, the isotope mass is not written out explicitly in the remaining sections. The 
final mixture ratio is determined in each experiment using infrared spectroscopy and previously 



determ ined CO, CO2 and H2O transmission band strengths with ~20% uncertainty dGerakines et al 



1995). The same band strengths are used to calculate absolute ice thicknesses, but in the thin ice 



experime nts the band str engths are first scaled to account for the longer absorption pathway in 
RAIRS dOberg et al.1120091) . which results in a ~50% uncertainty in the absolute ice thickness in 
these experiments. Based on the results of previous spectroscopic studies on mixed i ces, the H2O 



abundance is determined fr om the bending mode rather than the stretching feature dOberg et al 
2007llBouwman et al-booi . 



In each segregation experiment the ice mixture is deposited at a temperature T^ ep and then 
quickly heated (~5 K min" 1 ) to a chosen segregation experiment temperature r exp . The relative 
temperatures are controlled to a fraction of a degree, while the absolute temperatures have a ~2 K 
uncertainty. The segregation is measured by acquisition of infrared spectra during 2-4 hours; the 
spectral profiles of pure and mixed CO2 ice are distinctly different and can thus be used to measure 
the amount of segregated ice. In the thin ice experiment segregation is evaluated by using the CO 
and CO2 stretching bands and in the thick ices by combining information from the CO2 stretching 
and bending features. 

Table[T]lists the experiments carried out under ultra-high vacuum conditions. The ice-mixture 
thickness is varied between 8 and 37 ML to explore the thickness dependence in the ice-thickness 
regime accessibl e by RAIRS; the CO2 stretching fea ture is linear in absorbance for the first ~5- 



10 ML C0 2 ice (Te olis et al 



2007 



Oberg et al.ll2.009b . At higher coverage both the spectral shape 



and absorbance are affected by interference. This is only true for pure CO2, thus a 30 ML thick 
2:1 H20:C02 ice mixture is still accessible. Most H20:GC>2 experiments are carried out with a 
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Table 1. UHV ice-segregation experiments. 



Exp. 


C0 2 /CO 


H 2 : X 


Thick. (ML) 


Tdep (K) 


T exp (K) 


1 


C0 2 


10:1 


37 


19 


60 


2 


co 2 


3:1 


36 


19 


60 


3 


C0 2 


2:1 


8 


19 


53 


4 


co 2 


2:1 


8 


19 


58 


5 


co 2 


2:1 


11 


19 


50 


6 


co 2 


2:1 


11 


19 


55 


7 


co 2 


2:1 


11 


19 


60 


8 


co 2 


2:1 


13 


19 


50 


9 


co 2 


2:1 


18 


19 


55 


10 


co 2 


2:1 


18 


19 


60 


11 


co 2 


2:1 


20 


19 


53 


12 


co 2 


2:1 


27 


19 


56 


13 


co 2 


2:1 


30 


19 


50 


14 


co 2 


2:1 


17 


40 


55 


15 


co 2 


2:1 


21 


50 


55 


16 :l 


co 2 


2:1 


10 


19 


55 


L7 b 


co 2 


2:1 


9/75 


19 


53 


18 


co 2 


1:1 


13 


19 


50 


19 


co 2 


1:1 


13 


19 


55 


20 


CO 


2:1 


22 


19 


27 


21 


CO 


1:1 


10 


16 


27 


22 


CO 


1:1 


32 


19 


23 


23 


CO 


1:1 


31 


19 


25 


24 


CO 


1:1 


29 


19 


27 



"The ice was annealed at 60 K for less than a minute before cooled back to the segregation temperature. 

b The H 2 0:C0 2 ice mixture of 9 ML was deposited on top of 75 ML H 2 ice. The H 2 ice was deposited at 

100 K to achieve a compact ice substrate and then cooled to 19 K before depositing the ice mixtures. 

2:1 mixture, but four other mixing ratios are included as well (Exps. 1, 2, 18 and 19). The 2:1 
mixture was chosen as a standard experiment to facilitate the study of segregation dependencies 
over a range of conditions, which is not possible for more dilute mixtures at laboratory time scales. 
The deposition temperature is varied between 20 and 50 K and segregation is investigated between 
50 and 60 K. Two additional experiments (Exps. 16 and 17) explore the impact of fast thermal 
annealing and of using an ice substrate instead of depositing the ice directly on the gold surface. 
The ice substrate is constructed by depositing a thick layer of H2O ice onto the gold surface before 
depositing the ice mixture. The segregating ice mixture is thus isolated from the gold surface. Five 
experiments on H20:CO segregation (Exps. 20-24) at different temperatures and for different ice 
thicknesses and mixture ratios are included to test whether the H^OiCC^ segregation results can 
be generalized to other ice compositions. Table|2]lists the H20:CC>2 experiments carried out under 
high vacuum conditions, which are set up to explore temperature- and mixture-dependencies of ice 
segregation for thick ices between 100 and 300 ML and how this compares with the UHV thin ice 
experiments. 
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Table 2. HV ice-segregation experiments. 



Exp. 


H 2 : C0 2 


Thick. (ML) 


Tfep (K) 


i exp (K) 


I 


4:1 


430 


15 


60 


II 


4:1 


510 


15 


70 


III 


2:1 


140 


15 


40 


IV 


2:1 


160 


15 


45 


V 


2:1 


240 


15 


50 


VI 


2:1 


160 


15 


60 


VII 


2:1 


180 


15 


70 



3. Monte Carlo simulations 



Monte Carlo simulations are used to qualitatively investigate ice segregation through two different 
diffusion mechanis ms. The stochastic me tho d employed here builds o n the lattice-gas Monte Carlo 
technique used by 



Los et al 



(2006, 



2007) and 



Cuppen & Herbstl(l2007l) . In the simulations CO2 and 



H2O molecules are followed individually as a function of time, both position and environment, as 
they diffuse through the ice and desorb from it. The probability of a certain event in the simulation 
is governed by its energy barrier, which depends both on the kind of event and on the total energy 
difference before and after the event takes place. Exothermic events are therefore overall more 
likely to occur. Diffusion thus drives ice segregation when it is energetically favorable for molecules 
of the same kind to bind together rather than to bind with other types of molecules. The specific 
energy conditions for different types of diffusion is discussed further below. 

In the simulation, the ice structure is modeled as a regular lattice with dimensions 25 x 25 
(extended infinitely by periodic boundary conditions) x (/; + 1), where h is the ice thickness in 
monolayers. The lattice is initially filled randomly with the two molecular species investigated for 
segregation, such that the first (h - 1) layers are fully filled with molecules, the h layer is 75 % 
filled and the last h + 1 layer is 25 % filled. The partial fill of the top layers simulates the expected 
roughness of a surface. The ratio between the two molecular species occupying the sites is a given 
input parameter. 

The molecular distribution changes with time through two diffusion mechanisms, hopping of 
molecules into empty sites and swapping of molecules between sites, and through desorption into 
the gas phase. The rates for these three processes R^, R s and (where h is for hopping, s for 
swapping and d for desorption) are defined by 



R, 



v.xe 



-Ei/T a 



(1) 



where i=h, s or d, v; is the vibrat ional frequency of the molecule in its binding site, which is 



10 12 s 1 for physisorbed species jCuppen & Her bst 2007), T exp is the temperature at which the 



ice mixture is kept, and E\ is a barrier height in K for process i. 

All barrier heights depend on the difference in binding energy before and after an event. In 
a binary mixture of molecules A and B, the total binding energy E\, depends on the number and 
nature of (next)-nearest neighbors (n^ n , nf m , n^ nn , and surrounding the considered molecule 
and also on the specific molecule-to-molecule binding energies s£~ A , £j^~ B and s^~ B . The binding 
energy for next-nearest neighbors is defined to be i weaker than for nearest neighbors because of 
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the cubic geometry of the lattice. There is also an additional bulk energy term £t,uik, accounting for 
longer range effects. The total binding energy of a molecule A is thus described by 

F A — r-A-A^A _i_ "nnn\ , „B-A/„B , "nnn s 77 ,"->\ 
E b ~ B b ("nn + -g~) + S b ("nn + — ) + ^bulk- (2) 

The total hopping energy barrier for a molecule A into a neighboring site is defined by 



£,t = £h + ~A£ b A , (3) 



where e A is the barrier energy for a molecule to hop into a site where the total binding energy is the 
same before and after the event and A£ A is the difference in binding energy between the new and 
the old site. Hopping one site further has a barrier proportional to 2e A . The total swapping energy 
barrier is defined by 



+ -A£ A + -AE* (4) 

2 b 2 b> 



where e A ~ B is the barrier energy for two molecules A and B to swap into sites where the total 
binding energy is the same before and after the event. The desorption energy £ A equals the binding 
energy Zs A . 

With the rates and barriers thus defined, the Monte Carlo method is used to simulate which of 
these processes occurs for a specific molecule through the random walk method. The simulation 
starts with choosing time step At by picking a random number y between and 1 and calculating 
the event time Af event from 

-ln(y) 

Afevent " R h+ R s+ R A (5) 
for all molecules in the ice. The molecule with the smallest event time is selected and its event 
time defines the calculation time step. This molecule will then hop, swap or desorb at f culTen t + Af 
depending on the value of a second random number z between and + R s + /?d- The molecular 
configuration in the lattice is subsequently updated to start a new calculation cycle. 

In the simulation the energy barriers are chosen to represent a H20:C02 ice mixture (Table 
[3j. Since this is a qualitative model, the values are approximate. The binding energies are based 

H O— CO CO —CO 

on TPD experiments, except for e. 2 , which is assigned a lower value than s b - - to make 
CO2 segregation through hopping energetically favorable within the framework of the simulation. 
This may be counter intuitive since CO2 and H2O should form a hydrogen bond, which is stronger 
than the van-der-Waals interactions between two CO2 mol ecules. From TPD experiments there is 



H O— CO CO —CO 

however no evidence for e,~ 2 > s, ~ "in ices (Collings et al. 2004). This may be due to 



steric effects, i.e. CO2 molecules may be able to bind to more molecules in a pure CO2 ice compared 
to in a H20-rich environment, resulting in an overall stronger bond for a CO2 molecule that hops 
from a H2O dominated site to a CO2 dominated site despite the lower bond strength between 
individual CO2 molecules. This requires both more experimental and theoretical investigations. 
Until such data exist, segregation through hopping cannot be excluded and it is thus included in the 
simulations, though segregation through swapping alone is simulated as well. 

In contrast to CO2 hopping, swapping only requires that the average pure ice band strengths 
are greater than the average band strengths in the mixture to drive segregation. In other words 
for the swap of a CO2 molecule from an H2O dominated environment with a H2O molecule in a 
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it q Ti q CO —CO H O— CO 

CO2 dominated environment to be energe tically favorable e™ w " 2U + e^"^ > 2 x e™-" 2 . 
This is expected to hold, since H2O can form considerably stronger hydrogen bonds to other H2O 
molecules compared to with CO2. In addition hopping of H2O drives segregation under the condi- 
tion g^o-H^o > gfto-co, tn0U g n p rocess was j ess important in driving segregation at the tem- 
peratures and barriers chosen for the simulations. Thus diffusion in the ice will drive segregation 
under a range of different conditions with regards to relative energy barriers, and the dominating 
segregation mechanism may differ between different ice mixtures. 



The ho pping energies are multiples of the binding energies similarly to ICuppen & Herbst 



(2007) and Ga rrod et al.l d2008l) . while the swapping barrier is unknown and we chose to set it 
to a factor three higher than the CO2 hopping barrier. To investigate segregation through only hop- 
ping and only swapping, we also run simulations with only one of the mechanisms turned on. This 
is done computationally by increasing the barrier height by a factor of 10 for the excluded process 
compared to Table [3] 

Table 3. Defined binding, hopping and swapping energies of H2O and CO2. 



Barrier type a Energy / K 



H10-H10 

£ b~ 


1000 


C02-C02 

e b - - 


500 


HiO-COt 

V 


400 


H,0 
£ h~ 


2400 


co 2 

b h 


1200 


H,0-C0 2 


3600 


COi-HiO 


3600 



"b = binding, h = hopping, s = swapping. 



4. Results and analysis 

The results and analysis of H20:CC>2 segregation under UHV and HV conditions, H20:CO seg- 
regation under UHV conditions and Monte Carlo simulations are presented sequentially below. 
The CO2 UHV experiments form the central part of this section, while the other studies present 
supporting information about the segregation mechanism. 

4.1. UHV CO2 ice mixture experiments 

Figure [TJ shows the evolution of the CO2 stretching mode with time in a H20:CC>2 2:1 ice mixture 
at 55 K during four hours. Following deposition at 20 K, the ice is first heated to 50 K, where the 
ice mixture spectra is acquired, and then further to 55 K. At 55 K a new feature appears within a 
few minutes at 2310 cm -1 , which is attributed to pure CO2 ice (Fig. [2j. The ice segregates quickly 
during the first hour before the rate levels off; after four hours the segregation still proceeds slowly 
however. Similarly to this experiment, all investigated H20:C02 2:1 and 1:1 ice-mixture spectra 
change consistently with ice segregation when the ice mixtures are heated to 50-60 K and kept at 
the chosen segregation temperature for a minimum of two hours. 
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CO 

£> 

o 

CO 
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Fig. 1. Spectra illustrating the changes in the CO2 stretching feature as the ice segregation proceeds 
in a 18 ML thick 2:1 H20:C02 ice at 55 K (Exp. 9) over the course of 4 hours. The spectra are 
acquired every 5-15 minutes. The stretching band is shifted compared to normal CO2 because 
13 C0 2 is used. 

The early, fast segregation is inferred to be a surface process; the observed level of segregation 
does not depend on the initial ice thickness between 8 and 27 ML at early times for all ice temper- 
atures, and also at late times for experiments at 55 K or colder. Figure^ demonstrates this for 2: 1 
H20:CC>2 ice mixtures of three different thicknesses after three hours of segregation at 55 K, when 
the absolute amount of segregated ice is the same for all three ice mixtures within the experimental 
uncertainties. Segregation is somewhat thickness dependent at 60 K between 1 1 and 18 ML, which 
is indicative of a second bulk segregation process becoming efficient at this higher temperature 



Quantitatively, the amount of ice that is segregated at any time is calculated by simultaneously 
fitting a pure CO2 ice spectrum and an appropriate H20:CC>2 ice mixture spectrum - both acquired 
at 50 K, where segregation is still slow - to the segregating ice using an IDL script. During the 
automatic fit, the template spectra are allowed to shift with at most a few wavenumbers to obtain 
an optimal fit as examplified in Fig. [3] The fits are generally good except for a slight low-frequency 
mismatch, indicating that the structure of the mixed ice changes somewhat during heating. The re- 
sulting amount of segregated ice is plotted as a function of time in Fig.|4^ for three 13-18 ML thick 
2:1 H20:CC>2 ice mixtures at 50, 55 and 60 K, demonstrating that the segregation rate increases 
with temperature. The same trend is present for both thinner and thicker ices (not shown). 

The segregation rate also depends on the ice mixing ratio as shown in Fig. |4j? for 1:1 and 2:1 
H20:CC>2 mixtures at 55 K. The initial rate is ~3 times higher in the H20:CC>2 1:1 ice. This is also 
the case when comparing the 1 : 1 and 2: 1 mixtures at 50 K. The upper limit on the segregation rate 



(Fig. |2p). 
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Fig. 2. The spectra illustrate the level of segregation the ice reached in 2:1 H20:C02 ice mixtures 
of different thicknesses in a) after ~3 hours at 55-56 K and in b) after ~2 hours at 60 K. The ice 
thicknesses are 27 ML (black solid), 18 ML (blue dashed) and 11 ML (green dotted). The partly 
segregated ices are plotted together with a spectra of pure CO2 ice at 50 K (red dashed dotted line) 
in panel a). 



in the 10:1 mixture is a factor of 10 lower than the 2:1 rate. The segregation rate thus increases 
with a power 2-3 with the CO2 concentration in the ice mixture (where the concentration is with 
respect to the total H2O+CO2 ice amount). Furthermore the steady-state segregation level is 70% 
higher in the 1 : 1 ices compared to the 2: 1 mixtures, which is comparable to the CO2 concentration 
increase. This suggests that the same total amount of ice is available for surface segregation in all 
ice mixtures and that the final level of ice segregation only depends on the amount of CO2 in the 
top available ice layers. 
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Fig. 3. The spectra of an 18 ML thick H20:CC>2 ice with the initial composition of 2:1 (black) is 
fitted with a pure CO2 spectra (23 1 1 cm ) and a H^OtCCh 2: 1 mixture spectra at 50 K (red solid) 
after ~3 hours at a segregation temperature of 55 K. The combined fit is the red dashed line. 

In contrast to temperature and composition, the segregation rate is independent of the other 
experimental variables, including thermal annealing, deposition temperature below 50 K and the 
nature of the substrate (Figs. @t and d). Figure @t shows that an ice segregating at 55 K after 
a 1 min warm-up to 60 K behaves as a corresponding amorphous ice at 55 K, except for the 
first point, which can be attributed to segregation at 60 K. The same figure shows that increasing 
the deposition temperature from 18 to 40 K does not affect the segregation, while at 50 K the 
ice segregates during deposition (not shown). Any effect by the substrate on the segregation was 
investigated by comparing the segregation of two thin ices (8-9 ML) at 53 K after deposition 
directly on the gold substrate and on top of H2O ice; there is no measurable difference between the 
outcome of the two experiments. 

The initial segregation rate is measured quantitatively in each experiment by fitting an exponen- 
tial function A M (1 — e~W) to the first 40-240 minutes of the segregation, where A„ is a steady-state 
amount of segregated ice, t is the segregation time in seconds and £ seg is a temperature dependent 
rate coefficient describing the segregation process. Such a functional form is expected for a process 
that eventually reaches a steady-state and as demonstrated in Section 4, this is consistent with the 
simulations of segregation through surface diffusion. Fig. [4] shows that the equation fits the early 
segregation well but at later times, the experiments and exponential fits deviate, especially for the 
warmest ices (Fig. [4^, 60 K). This suggests two types of segregation mechanisms with different 
energy barriers, which are assigned to surface and bulk processes in Section 4. There are however 
too few points to fit both types of segregation in most experiments and thus we focus on the initial 
segregation rate alone. This is done by fitting the exponential growth function to early times, which 
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Fig. 4. The amount of segregated CO2 ice (in ML) as a function of time in H20:CC>2 ices. Panel a) 
shows the increasing segregation rate with temperature when comparing experiments at 50 (blue 
triangles), 55 (red diamonds) and 60 K (black stars) (exp. 8-10). Panel b) shows the segregation 
at 55 K in a 1:1 (red diamonds) and a 2:1 mixture (black stars) in two 13 ML thick ices (exp. 9 
and 20). Panel c) compares the segregation at 55 K for a quickly annealed ice (red diamonds) with 
amorphous ice deposited at the 'normal' 19 K and at 40 K (exp. 8, 14 and 16). Finally the curves in 
d) show the segregation at 53 K for 8-9 ML thick ices deposited directly on the gold substrate and 
on top of 75 ML compact H2O ice respectively (exp. 3 and 17). In each experiment, the segregation 
during the first 40-240 min is fitted with a function Aoo(l - e~ k * c "') (solid lines). 



is defined as 240 min at 50 K, 90 min at 52-56 K and 40 min above 56 K. At these times a single 
exponent fits the growth equally well as two exponents and the segregation is thus dominated by 
the fast process. To test the sensitivity of the fits to the fitting time intervals, these were varied by 
+50%. This does not affect the derivation of the segregation barrier, significantly, as shown further 
below. The results of these fits are shown for all experiments in Table |4] The rates are independent 
of ice thickness for 8-20 ML thick ices, while the thickest ice mixtures segregate somewhat slower 
compared to thinner ices kept at the same temperature. This decrease is barely significant, however, 
and it is probably due to RAIRS effects, which become important at these thicknesses as discussed 
in Section 2, rather than to real differences in the segregation behavior above 20 ML. 

Assuming that the segregation-rate coefficients can be described by the Arrhenius equation, 
k — T x e~ £ » g /7 where T is the pre-exponential factor, E seg the segregation barrier in K and T the 
ice temperature, the segregation energy can be extracted by plotting ln(fc) versus l/T and fitting a 
straight line to the points; the slope of the line is E seg . This is shown in Fig.[5]for the rate coefficients 
belonging to the 8-13 ML experiments, resulting in a segregation barrier of 1080 + 190 K for the 
2:1 and 1050 + 220 K for the 1:1 H20:CC>2 ice mixtures, where the uncertainties include both 
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Fig. 5. The natural logarithm of rate coefficients for segregation in thin ices (8-13 ML) is plotted 
versus the inverse of the segregation temperature for H20:C02 2: 1 (black crosses) and 1 : 1 mixtures 
(red diamonds). The black line is fitted to the 2:1 experiments. 
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experimental and fit errors. If the time intervals chosen to fit the first exponential is allowed to vary 
by +50% from the optimal value, the derived barrier is instead 1200+250 K for the 2:1 mixtures. 
This probably overestimates the uncertainties. Still it is comforting that the derived value is well 
within the uncertainty of the segregation barrier derived from the optimal fitting time intervals, and 
the presented barrier for surface segregation is thus robust. The barriers are also not significantly 
different for the two mixtures, even though the rates are higher for the 1:1 ices; the pre-exponential 
factors amount to 2 x 10 (5±I) for the 2:1 mixtures and 6 x 10 (5±I) for the 1:1 mixtures. The pre- 
exponential factor can be understood as a diffusion frequency multiplied by the fraction of diffusion 
events that results in segregation. Its increase between the two mixing ratios is thus consistent with 
the increase in the number of CO2 molecules available in the top layers. 

Table |4] shows that the steady-state amount of segregated ice from the fast segregation mecha- 
nism is consistently between 1 and 2 ML in all experiments, except for the 1 : 1 mixtures as discussed 
above. The steady-state amount increases slightly between the 50-55 K and the 60 K experiments. 
This increase is however low enough that an average value over all ice temperatures and thick- 
nesses can be used when quantifying the segregation process. Combining the information in this 
section, the amount of segregated ice re™ in ML after a time t in seconds at an ice temperature T 
in K is well described by 



where x is the CO2 concentration in the ice with respect to the total CC^FbO ice amount (in 
astrophysical settings this will be dominated by the H2O ice fraction). At 60 K the second, bulk 
segregation rate is an order of magnitude lower than the initial surface segregation rate. There are 
too few measurements to quantify the bulk segregation further, but this difference suggests that the 
barrier to segregate in the bulk is more than a factor of two larger than for the surface segregation 
evaluated here. In Section 4.4 this bulk segregation is simulated through molecular pair swapping 
with a three times higher energy barrier compared to CO2 surface hopping. 

4.2. HV CO2 ice mixture experiments 

Transmission spectroscopy in the HV set-up allows for the investigation of segregation in an or- 
der of magnitude thicker ice than is possible in CRYOPAD. Figure [6] shows the segregation of 
H20:C02 2:1 mixtures through spectroscopy of both the CO2 bending and stretching modes. Some 
ice rearrangement is also visible in the 4:1 mixture but there is no clear indication of segrega- 
tion into two different ice phases. Thus, similarly to the thin ice experiments, the segregation rate 
depends on the ice mixture ratios. 

The amount of segregated ice at each time step is derived using the same procedure for the 
thick ices as described above for the thin ices. Figure [7] shows the fraction (not the number of ML) 
of the ice mixture that is segregated as a function of time at segregation temperatures between 40 
and 70 K. Above 50 K the segregation of the >140 ML ices reaches 50%, from which we infer, in 
agreement with previous studies, that segregation is a bulk process for thick ices deposited under 
high-vacuum conditions. It also shows however that a large fraction of the ice is protected from 
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Fig. 6. CO2 bending (~635 cm 1 ) and stretching (~2270 cirT 1 ) spectral features at the onset (black 
solid) and after 2 hours of heating at the segregation temperature (red dashed) for three 160-510 
ML thick H20:CC>2 ice mixture experiments. 
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Fig. 7. The segregated ice fraction, with respect to the total CO2 ice amount, as a function of time 
at 40 (black stars), 50 (blue triangles), 60 (red diamonds) and 70 K (green circles) for 140-240 ML 
thick H20:CC>2 2: 1 ices. The solid lines show the exponential fits from which the rate coefficient is 
derived for each experiment. 
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segregation at all ice temperatures, though this fraction is smaller than in the thin ice experiments 
that are dominated by surface diffusion. The segregated amount of ice as a function of time is fit 
to an exponential function in the low-temperature experiments - above 45 K, segregation is too 
fast for a reliable rate determination, i.e. the final segregated amount is reached within a couple 
of measurements. In contrast to thin ices the segregation can be fitted with a single exponential 
function for all temperatures and the determined steady-state amount of segregated ice depends 
on the ice temperature; the segregated fraction varies between 20% and 55% at the investigated 
temperatures. The resulting bulk segregation rates at 40 and 45 K are (2 + 6)x 10~ 4 and (5 + 3)x 10~ 3 
s _1 , respectively. Since the rates describe the segregation of a fraction of the ice rather than a 
number of monolayers they cannot be directly compared to the thin ice surface segregation rates - 
it is, however, clear that steady state is reached in a shorter time in the thick ice experiments than 
in the thin ices. The derived barrier from the 40 and 45 K experiments is 1270 + 560 K, which is 
consistent with the thin ice segregation barrier of 1080+190 K. The thick ice segregation, while 
more efficient, thus does not necessarily proceed through a separate mechanism. 



4.3. UHV CO ice mixture experiments 



Because of the high volatility of CO, H20:CO segregation is investigated between 23 and 27 K. 
This is below the temperature range for the H2O ice transitio n between high and low den sity amor- 



phous phases, which occurs gradually between 38 and 68 K (IJenniskens & Blake 



1994). The tem- 



peratures investigated for H20:CO and H20:CC>2 segregation are similar in terms of % of the CO 
and CO2 desorption temperature; CO desorbs at ~30 K and CO2 at ~70 K from a H2O ice sub- 
strate. Hence, if segregation occurs through a simple diffusion process with barriers proportional to 
the binding energy in CO and CO2 ice mixtures, this is the predicted range of temperatures where 
the H20:CO ice should segregate. Figure [8] shows the H20:CO ice evolution when quickly heated 
from 16 to 25 K and then maintained at 25 K for four hours. The new band appearing around 2094 
cm 1 at 25 K has the same width and a comparable position to a pure CO ice feature and this is 
interpreted as resulting from ice segregation. Segregation thus occurs and it occurs only in a part 
of the ice, similarly to the CO2 mixtures. 

The segregation is quantified by fitting three Gaussian's to the spectral feature, two belonging 
to CO ice in a H2O mixture and one corresponding to the pure CO ice band (Fig. [9]). The two 
mixed-ice Gaussians have bee n observed previously in several studies and are attributed to two 
different CO-H2O interactions dCoilings et al.ll2003l) . The Gaussian parameters are free variables 
during the fit except for the FWH M of the pure CO Gauss ian, which is constrained to be within 0.5 
cirT 1 of the pure CO ice value in iBouwman et alj d2007l) . The resulting Gaussian parameters for 



all H20:CO experiments are reported in Table [5] The wid ths of Gaussian s 1 and 3 are similar to 

J2OO7J). Figure |9] shows 



Bouwman et al 



those used to fit the same H20:CO ice-mixture feature by ] 
that the spectral fit cannot exclude a few % of ice segregation at deposition. This does not affect the 
results since the reported level of segregation during the experiment is calculated from the increase 
of Gaussian 2 compared to deposition. 

Using the Gaussian fits, Fig. [10] shows the time evolution of the amount of segregated ice 
in the different CO ice mixtures, except for the 23 K ice where no segregation is observed. The 
segregation time series have the same shape as the CO2 segregation curves, with a fast initial 
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Fig. 8. Ice segregation in a 10 ML H20:CO mixture observed through a) the changes in the CO 
stretching feature with time (from blue to orange), plotted together with b) the original ice mixture 
spectra (black solid) and c) a pure CO ice spectra (black dashed). The spectra are off-set compared 
to the normal isotope band positions, because I3 CO is used. 
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Fig. 9. The resulting fit of three Gaussians to a) a partially segregated 10 ML H20:CO ice at 25 K 
and b) the original ice mixture at 16 K. The green dotted Gaussians correspond to mixed ice fea- 
tures, the blue dashed-dotted to pure CO ice and the combined fit is plotted in red dashed lines. 
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Table 5. Parameters for the three Gaussians used to fit the CO feature. 
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Fig. 10. The amount of segregated ice as a function of time in four H20:CO mixtures: a) 1 : 1 , 27 K, 
29 ML, b) 2:1, 27 K, 22 ML, c) 1:1,25 K, 31 ML and d) 1:1, 27 K, 10 ML. 

segregation resulting in <2 ML segregated ice, followed by a slower segregation process at the 
highest temperature. The ice composition dependence seems weaker, but is still present, compared 
to the CO2 case. The segregation at 27 K is thickness dependent, which suggest that bulk diffusion 
has already become important at this temperature, though the initial rate should still be dominated 
by surface processes. Comparing the early-time rates of the 1:1 mixtures at 25 and 27 K results 
in an approximate barrier for the initial H20:CO segregation of 300 + 100 K, using the same 
fitting technique as for the CO2 mixtures. This is significantly lower than the H20:CC>2 segregation 
barrier, as would be expected for the more volatile CO ice. 

4.4. Monte Carlo simulations 

This section presents the results from six different ice segregation simulations, investigating the 
importance of swapping and hopping, and the effect of temperature on the segregation behavior. 
Figure QT| shows the simulated H20:C02 ice segregation in a 5 ML thick 1 : 1 mixture using the 
Monte Carlo technique and inputs described in Section 3. The segregated ice fraction is quan- 
tified through an increase in the number of nearest neighbors of the same species, which is not 
directly comparable to the spectroscopic segregation observed in the experiments, though the two 
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are linked. The simulated amount of segregated ice is thus presented in arbitrary units. With the 
specified energy barriers, segregation at 55 K occurs through a combination of hopping and swap- 
ping; removing one of the mechanisms decreases the total segregated fraction (Fig.fTTT). 

The segregation through hopping is well defined by an exponential function, i.e. it is a single 
process approaching steady-state, while segregation through swapping, and segregation through 
swapping + hopping are both not. Instead segregation through swapping is initially fast and then 
continues to proceed slowly until the end of the simulation, long after the segregation through hop- 
ping has reached steady-state. The initial fast swapping is due to swapping with surface molecules, 
which explains that segregation through hopping and swapping is not additive. It also explains 
why segregation through only swapping cannot be fitted by a single exponential function, since 
it proceeds through two different processes, surface and bulk swapping. Of the three curves the 
swapping+hopping one is qualitatively most similar to the 55 and 60 K CO2 UHV experiments 
and hopping alone cannot explain the non-exponential behavior of the experimental segregation 
curves during late times. It is however important that swapping alone may be able to reproduce 
the segregation behavior observed in the experiments and this implies that the CO2-CO2 binding 
energy is not required to be higher than the CO2-H2O one for segregation to occur through a com- 
bination of surface and bulk diffusion. 

The segregation rate increases with segregation temperature in the simulation when both hop- 
ping and swapping are included, as is observed in the experiments. Figure [12] also shows that 
the curves at higher temperatures, i.e. 55 and 60 K, diverges more from an exponential growth 
compared to the 50 K simulation, suggesting an increasing importance of bulk-swapping with tem- 
perature. This is also consistent with the experimental results (Fig.|4^). The segregated fraction in 
the simulations increases more clearly with temperature than in the experiments, suggesting that 
the bulk process, swapping of bulk molecules, may be somewhat too efficient in the simulations. 



5. Discussion 

Following a brief comparison with previous experiments, this section combines the results from 
(U)HV H20:C02 and H20:CO segregation experiments and Monte Carlo simulations to constrain 
the segregation mechanics and the meaning of the empirically determined segregation rates and 
barriers. This information is then used to discuss the astrophysical implications of quantifying the 
segregation of binary ice mixtures. 



5. 1 . Comparison with previous experiments 



Previous studi es have focused on segregation in thick ices (>100 ML) under HV conditions. Of 
these studies, iHodvss et alj (|2008) carried out the mos t similar expe r iment s to our HV ones and 

(2008) reported an increase 



Hodvss et al. 



thus we only compare these two sets of experiments, 
in the segregated ice fraction and in the segregation rate between 55 and 70 K for a H20:C02 4: 1 
mixture, which is comparable to what we observe between 40 and 60 K for a 2: 1 mixture. The lower 
temperatures required in our experiment to achieve the same result are co nsistent with the obse rved 
dependence on composition for segregation in both studies. Similarly to 



Hodv ss et al. 



(2008), we 



do not observe any dependence on the ice deposition temperature in our thin ice experiments, as 
long as the deposition temperature is lower than the lowest temperature for which segregation is 
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Fig. 11. Simulations of ice segregation in a P^OrCCh 1:1 mixture at 55 K due to a) swapping and 
hopping, b) only hopping, and c) only swapping. The dashed lines show exponential fits to the first 
30 min. 



observed. Qualitatively, the studies also agree on the segregation dependence on composition. A 
quantitative comparison is not possible since such data have not been reported previously. 

5.2. Segregation mechanisms 

The possible segregation mechanisms for thin and thick ice mixtures are summarized in Fig. Q~3] 
Thin binary H20:C02 ice mixtures, 8-37 ML, segregate differently compared to thick, >100 ML, 
ices. This implies either a thickness dependent segregation mechanism or a different ice structure 
above 100 ML compared to <40 ML. Segregation in thin ices is qualitatively reproduced by Monte 
Carlo simulations including surface hopping, surface swapping and the slower molecular swapping 
throughout the bulk of the ice. If thicker ices are more porous because of cracking, hopping and 
surface swapping may be pos sible throughout the ice resulting in the observed fast bulk segregation. 



Previous studies, including iHodvss et alj (|2008), instead suggest that segregation is the product 



of a H2O phase transition around the segregation temperature. Phase transitions may require a 
certain ice thickness to become effective and thus may only be important for thicker ices; the 
H20:CO results show that phase transitions are not responsible for segregation in thin ices. The 
similar segregation barrier in thick ices and thin ices suggests that the mechanism is the same at 
all thicknesses and that the rate differences have a structural origin. The shorter times required 
to reach steady-state in thick compared to thin ices indicates, however, the opposite. More thick 
ice experiments containing molecules with different segregation temperatures are thus required to 
differentiate between these two segregation scenarios. In either case, these different segregation 
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Fig. 12. Simulations of ice segregation in a H20:CC>2 1:1 mixture at 50, 55 and 60 K due to a 
combination of swapping and hopping. The dashed lines show exponential fits to the first 30 min. 



processes in thin and thick ices show that ice dynamics, including ice chemistry, is not necessarily 
thickness independent as has been previously assumed in the astrochemical literature. 

The segregation growth curves of thin H20:CC>2 and H20:CO ices have the same shape with 
an initial exponential growth up to 1-2 ML, followed by a slower growth at longer time scales. 
The segregation temperatures and inferred barriers differ by a factor 2-3. T he desorption/bindin g 



energies of pure CO2 and CO ices are similarly a factor of 2-3 different dCollings et al 



2004). 



The derived segregation barriers for the two types of mixtures are thus consistent with segregation 
through molecular hopping and swapping, with hopping and swapping barriers proportional to the 
desorption/binding energy of the most volatile mixture component. The similar ratio between the 
surface segregation barrier and desorption energies in the two types of ice mixtures, suggests that 
ice surface dynamics can be approximately parameterized based on desorption energies of pure 
ices alone. 

Physically, it is not obvious from the simulations whether this initial surface segregation is dom- 
inated by surface hopping or surface swapping or whether it is even the same in both the FbOrCC^ 
and H20:CO mixtures. Swapping is certainly required to explain the long-time growth of the seg- 
regated ice fraction. Depending on the relative barriers for hopping and swapping and the different 
binding energies in the ice, it may also be responsible for at least a part of the surface segregation as 
shown in Fig.QT] The similar behavior of the H20:CC»2 H20:CO ices actually suggests that surface 
swapping rather than hopping drives surface segregation in both ices, since there i s experimental 



evide nce for a higher H2O-CO binding energy compared to the CO-CO one in ices (ICollings et al 



2003) and surface hopping can only drive segregation if the hopping molecules forms stronger 



bonds with its own kind. Quantitative H20:C02 desorption data is however required to exclude 
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Fig. 13. The segregation mechanism in a) ices below 40 ML together with the two proposed segre- 
gation mechanisms for thick ices (>100 ML), through b) internal surface segregation and c) phase 
transitions. Black indicates segregated CO2 ice and stripes the H20:CC>2 ice mixture. Each step 
is marked with the approximate temperature at which the mechanism becomes important at the 
time-scales of low-mass star formation, as discussed in Section 5.3. 



surface segregation through CO2 hopping. A more quantitative simulation study of ice segregation 
is currently being pursued, which will address this question (Cuppen et al. in prep.). Until then, 
Eqs. 6 and 7 describe a non-specified surface segregation, while bulk segregation is approximately 
an order of magnitude slower in thin ices. From the experimental studies and the simulations, seg- 
regation is expected to occur in all ices, through swapping, where the binding energies of A-A 
+ B-B are greater than 2xA-B. The segregation barrier is expected to be mixture dependent, but 
related to the pure ice binding energies. 



5.3. Astrophysical implications 

Nature is 'non-cooperative' in that the typical astrophysical ice thicknesses of <100 ML are not 
guaranteed to be in the thin ice regime where hopping and swapping alone are responsible for ice 
segregation. Yet, since 100 ML is the upper limit assuming spherical grains, most astrophysical 
H20:CC>2 ice mixtures are probably closer to the thin ice regime than the thick ice one. Therefore 
this section will discuss ice segregation during star formation under the assumption that this can be 
described by the same equations as thin ice surface segregation. 
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During star formation ices are heated in the collapsing envelope to a certain temperature, de- 
pending on their distance from the protostar. The amount of time they spend at these elevated 
temperatures before desorbing depends on the infall rate that transports ices from the outer to the 
inner parts of the protostellar envelope. This time scale T cr ; t for ice heating in the envelope during 
low-mass star formation thus depends on several different parameters of which the stellar lumi- 
nosity is the most important one and will vary between different objects. Building on Eq. 7, the 
characteristic temperature r seg at which segregation is important can be calculated as a function of 
this heating time scale r cl -i t = k~ x for a known ice composition 



T 

1 seg 



1080 ± 190 



log^XlO^xfgfjf'^XTcrit) 



(8) 



where x is the CO2 concentration in the HyO i CO? ice phase, which is s et to 0.16 from the 5:1 ice 



mixture observed bv lPontoppidan et al 



(2008). 



Pontoppidan et al 



(2008) also calculated that under 



the assumption of free-fall accretion, a dust grain is heated to 25-50 K by a s olar-mass protosta r 
for ~4000 years, while in a recent infall model it is at least a factor of 5 longer ( iVisser et al I2OO9I 
Taking the lower limit 4000 yrs time scale for segregation, r seg is 30+5 K. This is low compared 
to previous estimates, of 50-100 K, and comparable to the temperature for CO distillation from 
CO:CC>2 ice mixtures, which is a competitive mechanism to form pure CO2 ice. With this segrega- 
tion temperature, segregation of H20:C02 ice is no longer excluded as a source for pure CO2 ice 
around low-mass protostars. This temperature is only for surface segregation however. Therefore if 
more than 2 ML of pure CO2 is present t owards a low-mass protostar, distillation is still required. 
Observations show that in the sample of iPontoppidan et alj (12008b . the average ice abundance of 
pure CO2 with respect to H2O ice is 2.4+2.4%. Similar values are also found towards the low-mass 
and high-mass sources in the sample when calculated separately for the two types of objects. With 
a total ice thickness of less than 100 ML, the amount of segregated CO2 towards most sources is 
thus consistent with surface segregation as long as the observations are dominated by ice at > 30 K. 

To understand the origin of pure CO2 ice towards individual objects requires estimates of the 
infall rate, the stellar luminosity and the origin of the ice absorption features in the envelope, which 
is outside of the scope of this study. Such a model could however use pure CO2 observations to 
constrain the thermal history of the protostellar envelope. While the characteristic temperatures of 
CO2 segregation in a H20:C02 ice mixture and CO distillation from a C02:CO ice are too close 
to differentiate their relative contributions towards most objects, their similarity is actually an asset 
when using pure CO2 ice as a thermal probe. Both types of ice processing are irreversible; hence 
once segregated or distilled, the CO2 ice will remain pure. Modeling of the amount of pure CO2 ice 
in a protostellar envelope may then reveal the maximum size of the envelope that has been heated to 
25-30 K. This inferred heating of the envelope in the past can be compared to the current luminosity 
of the protostar. A statistic on the discrepancy between the two would constrain the variability of 
protostars during an otherwise observationally inaccessible period of the stellar evolution. 

This picture may be complicated by the presence of strong UV fields around some protostars. 
UV radiation can produce CO2 from CO and H2O ices through photochemistry, photodesorb pure 
CO2 ice, and photodissociate CO2 into CO and thus transform the pure CO2 ice into a CQ2:CO 

; lLoeffler etalj|2005 : 



mixture (Watanabe & Kouchi 2002; 



Ober g et al 



2009; 



Gerakines et al. 1996). 



Of these processes, CO2 formation through photochemistry should not affect the amount of pure 
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CO2 ice, since CO2 then forms mixed with H2O and CO ices. UV-induced destruction of some of 
the produced CO2 ice may be important, however. Assessing its impact for a certain object requires 
an individual protostellar envelope model that includes UV-induced ice processes. Without such a 
model, pure CO2 ice observations still provides a lower limit to the amount of heat experienced by 
the ices in a protostellar envelope. 

6. Conclusions 

Quantifying ice dynamics is possible through a combination of a relatively large set of laboratory 
experiments, spanning the available parameter space, rate-equation modeling and Monte Carlo sim- 
ulations. This process is more time consuming than the qualitative studies we have been building 
on, but has the critical advantage of providing actual rates, which can be incorporated into astro- 
chemical models. These rates can then be used to test hypotheses about ice dynamics in space. 
Quantified ice processes can also be used as quantitative probes, and thus provide powerful tools 
to investigate processes during star formation. 

For the specific dynamical processes investigated in this study our main conclusions are: 

1. Thin (8-37 ML) H20:C02 and F^OrCO ice mixtures segregate through surface processes fol- 
lowed by an order of magnitude slower bulk diffusion. 

2. The thin ice segregation process is qualitatively reproduced by Monte Carlo simulations where 
ices segregate through a combination of molecule hopping and swapping. Segregation is ex- 
pected to be a general feature of ice mixtures where binding energies of all the mixture con- 
stituents are higher on average in the segregated ices compared to in the ice mixtures. 

3. The surface segregation barrier is 1080 + 190 K for the H20:C02 ice mixtures and 300 ± 100 K 
for H20:CO segregation. 

4. H20:C02 surface segregation is ice-mixture ratio dependent, quantified as k seg = 2 x 10(5*1) X 
e -i080±i90/r x ^_*_^ 25±1 ^ w here fc seg is the rate at which the top few monolayers segregate, T the 
ice temperature and x the CO2 concentration in the ice with respect to the total ice amount. 

5. Segregation in thick (> 100 ML) ices involves the bulk of the ice, which can be explained by 
ice cracking or H2O phase transitions or both. Ice dynamics thus depends on the thickness of 
the investigated ice. 

6. During low-mass star formation the surface segregation temperature of a typical thin HbOiCCh 
5:1 ice mixture is reduced to ~30 K, because of the longer time scales. 
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